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Abstract 



We consider the problem of solving dual monotone inclusions involving sums of composite 
parallel-sum type operators. A feature of this work is to exploit explicitly the cocoercivity of 
some of the operators appearing in the model. Several splitting algorithms recently proposed 
in the literature are recovered as special cases. 
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1 Introduction 

Monotone operator splitting methods have found many applications in applied mathematics, e.g., 
evolution inclusions [2], partial differential equations [1, 20, 23], mechanics [21], variational inequal- 
ities [6, 19], Nash equilibria [8], and various optimization problems [7, 9, 10, 14, 15, 17, 25, 29]. 
In such formulations, cocoercivity often plays a central role; see for instance [2, 6, 11, 13, 19, 20, 
21, 23, 28, 29, 30]. Recall that an operator C: Ti ^ Ti is cocoercive with constant (3 G ]0, +oo[ if 
its inverse is /3-strongly monotone, that is. 
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In this paper, we revisit a general primal-dual splitting framework proposed in [16] in the presence 
Lipschitzian operators in the context of cocoercive operators. This will lead to a new type of 
splitting technique and provide a unifying framework for some algorithms recently proposed in 
the literature. The problem under investigation is the following, where the parallel sum operation 
is denoted by □ (see (2.4)). 

Problem 1.1 Let be a real Hilbert space, let z £ Ti, let m be a strictly positive integer, let 
{coi)i<i<m be real numbers in ]0, 1] such that Ya^i = 1, let ^: "H — 2^ be maximally monotone, 
and let C : 7i ^ Ti he ^-cocoercive for some /i G ]0, +oo[. For every z G {1, . . . , m}, let Qi be a real 
Hilbert space, let rj G Qi, let 3^: Gi ^ 2^' be maximally monotone, let Di: Qi ^ 2^' he maximally 
monotone and z^'j-strongly monotone for some Vi G ]0,+oo[, and suppose that Li : Ti ^ Gi is a 
nonzero bounded linear operator. The problem is to solve the primal inclusion 

m 

find X Gn such that z e Ax + '^uJiL*[{Bi □ Di){LiX - r^)) + Cx, (1.2) 

1=1 

together with the dual inclusion 



find f 1 G ^1, . . . , £ Qra such that (3 x G 7^) 



z - YT=i ^i^i^i ^Ax + Cx 
^ (Vi G {1, . . . , m]) Vi G {Bi □ Di){LiX - r^). 

(1.3) 

We denote by V and V the sets of solutions to (1.2) and (1.3), respectively. 



In the case when {D- ^)i<i<m and C are general monotone Lipschitzian operators. Problem 1.1 
was investigated in [16]. Here are a couple of special cases of Problem 1.1. 



Example 1.2 In Problem 1.1, set z = and 

iO) and Dr. V ^ I 

if v 7^ 0. 



(Vi G {1,... ,m}) Bi:v^{{)} and Di-.v^T' ^ °' (1.4) 



The primal inclusion (1.2) reduces to 

find X & % such that Q ^ Ax + Cx. (1.5) 
This problem is studied in [2, 11, 13, 17, 23, 28, 29]. 

Example 1.3 Suppose that in Problem 1.1 the operators (Di)i<i<m are as in (1.4), and that 

A:x^{<d} and C:xh^O. (1.6) 
Then we obtain the primal-dual pair 

m 

find X G Ti such that z G ^iL* [Bi{LiX — r^)) , (1-7) 

i=l 

and 

(3 X G n){yi G {1, . . . , m}) Vi G Bi{LiX - r^). 
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This framework is considered in [7], where further special cases will be found. In particular, it 
contains the classical Fenchel-Rockafellar [27] and Mosco [24] duality settings, as well as that of 
[3]. 

The paper is organized as follows. Section 2 is devoted to notation and background. In 
Section 3, we present our algorithm, prove its convergence, and compare it to existing work. 
Applications to minimization problems are provided in Section 4, where further connections with 
the state-of-the-art are made. 



2 Notation and background 

We recall some notation and background from convex analysis and monotone operator theory (see 
[6] for a detailed account). 

Throughout, Ti^ G, and {Qi)i<i<m are real Hilbert spaces. The scalars product and the associ- 
ated norms of both H and G are denoted respectively by (• | •) and || • ||. For every i G {1, . . . , m}, 
the scalar product and associated norm of Gi are denoted respectively by (• | ■)g. and || • \\g.. We 
denote by B{H,Q) the space of all bounded linear operators from % to Q. The symbols ^ and 
— )• denote respectively weak and strong convergence. Let A: T-L ^ 2^ he a set- valued opera- 
tor. The domain and the graph of A are respectively defined by dom^ = |x S ^ | Ax ^ 0} 
and graA = {{x,u) G V. x V. \ u G Ax^ We denote by zev A = {x G | G Ax} the set of 
zeros of A, and by ran^ = {u £ H \ (3 x G H) u G ^x} the range of A. The inverse of A is 
A~^ : H ^ 2'^ : u {x e H \ u e Ax]. The resolvent of A is 

Ja = {U+A)-\ (2.1) 

where Id denotes the identity operator on Ti. Moreover, A is monotone if 

{y{x,y) xn) {y{u,v) e Ax X Ay) {x-y\u-v)>0, (2.2) 

and maximally monotone if it is monotone and there exists no monotone operator B: % ^ 2^ 
such that graB properly contains gra^. We say that A is uniformly monotone at x G domA if 
there exists an increasing function (p: [0, -|-oo[ — t- [0, -|-oo] vanishing only at such that 

(Vu G Ax)(V(y,i;) G gra^) {x - y \ u - v) > (j){\\x - y\\). (2.3) 

If ^ — aid is monotone for some a G ]0, +oo[, then A is said to be a-strongly monotone. The 
parallel sum of two set-valued operators A and B from l-L to 2^ is 

AUB = {A-^ +B-^)-^. (2.4) 

The class of all lower semicontinuous convex functions /: ^ — t- ]— oo,-|-oo] such that dom/ = 
{x G "^^ I /(x) < +00} / is denoted by ro(^). Now, let / G Tq{T-L). The conjugate of / is 
the function /* G Tq{H) defined by /* : u 1— )• sup^g-^((x | u) — /(x)), and the subdifferential of 
/ G TQi(H) is the maximally monotone operator 

df:n^2^: {uen\{yy (^U) {y - x \ u) + f{x) < f{y)] (2.5) 
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with inverse given by 



(df)-' = dr. 



Moreover, the proximity operator of / is 



1, 



proxj- : 7^ — )• 7^ : X I— )• argmin /(y) H — Hx — . 
We have 

Jdf = proxj . 

The infimal convolution of two functions / and g from % to ]— oo, +cxd] is 

f U g: H ^ ]— CO, +oo] : x i— )• inf (/(x) + g{x — y)). 

yen 



(2.6) 

(2.7) 
(2.8) 
(2.9) 



Finally, let 5 be a convex subset of %. The relative interior of S, i.e., the set of points x £ S such 
that the cone generated by x + S is a vector subspace oiTi, is denoted by ri S. 



3 Algorithm and convergence 



Our main result is the following theorem, in which we introduce our splitting algorithm and prove 
its convergence. 



Theorem 3.1 In Problem 1.1, suppose that 



G ran (a + Y,uj,L*{{Bi □ Di){U ■ -n)) + c\ 
^ i=i ^ 



(3.1) 



Let T and ((Ji)i<j<m be strictly positive numbers such that 



2pmin{/i, i/i, . . . ,Um} > 1, where p = min |r \ cr^ \ . . . ,o-„^| ~ r ^ (JiWi||LiP j . (3.2) 

Let e G ]0, 1[, let (\n)n&i be a sequence in [e,l], let xq G V., let (ai_„.)nGN cind (a2,n)neN be 
absolutely summable sequences in %. For every z G {1, . . . ,m}, let Vi^ G Qi, and let {bi^n)nm o-nd 
(cj,n)neN absolutely summable sequences in Qi. Let {xn)n<m and (f i,n> ■ ■ ■ , Vm,n)nm be sequences 
generated by the following routine 



(Vn G N) 



Pn = JtaIxu - t(^T^1Li ^iL*Vi^n + CXn + ai,„ - Z^j + a2,n 
Vn — 

for z = 1, . . . , m 

= '^(TiB^^ {vi,n + O'ii^Liyn — ^Vi^n ~ Ci^n — ^i^^ ~^ bi,n 
Vi,n+1 = Vi^n + K{qi,n " ««,«)• 



(3.3) 



Then the following hold for some x gV and (f i, . . . , Vm) G V. 
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(i) Xn^ X and (vi^n, ■ ■ ■ ,Vm,n) ^ (^^1, • • -^Vm)- 

(ii) Suppose that C is uniformly monotone at x. Then 

(iii) Suppose that DJ^ is uniformly monotone at vj for some j G {1, • • • ,m}. Then vj^n — ^ vj. 

Proof. We define Q as the real Hilbert space obtained by endowing the Cartesian product Qi x 
... X Qrn with the scalar product and the associated norm respectively defined by 



{■ \ ■)g : {v,w) ^'^uji{vi\ Wi)g_ and 



\ i=i 



where v = (vi, . . . ,Vm) and w = (wi, . . . , Wm) denote generic elements in Q. Next, we let K be 
the Hilbert direct sum 

K = neG. (3.5) 

Thus, the scalar product and the norm of /C are respectively defined by 



{■ \ ■)k ■ {(^^'")^iy^'^)) ^ \ y) + \ '^)g and W ■ Wfc- {x,v) ^ ^^\\xf + \\v\\l. (3.6) 
Let us set 

M : K^2^ 

ix,Vi,.. . ,Vra) ^ { - z + Ax) X {ri + B^'^Vi) X ...X {rra + B'^Vm)- (3.7) 

Since the operators A and (-Bi)i<i<m are maximally monotone, M is maximally monotone [6, 
Propositions 20.22 and 20.23]. We also introduce 

S-.K^K (3.8) 

{x,Vi, ... ,Vm)'-^ i ^UJiL*Vi, -Lix, -LmXj. (3.9) 



Note that 5^ is linear, bounded, and skew (i.e, S* = —S). Hence, S is maximally monotone [6, Ex- 
ample 20.30]. Moreover, since domS* = /C, M + S is maximally monotone [6, Corollary 24.24(i)]. 
Since, for every i G {1, • • • ,m}, Di is z/j-strongly monotone, D^^ is z^j -co coercive. Let us prove 
that 

{x,vi,..., Vm) ^ {Cx, D^'^vi, . . . , D:;:^Vm) (3.10) 

is /3-cocoercive with 

/3 = min{/i, i/i, . . . (3.11) 
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For every {x,vi, . . . , Vm) and every {y,wi, . . . , Wm) in /C, we have 

((x,t>i, ...,Vm)- {y,Wi, . . .,Wm) \ Q{x,Vi, ...,Vm)- Q{y,Wi, . . ■,Wm))jc 

m 

= {x - y \ Cx - Cy) + ^ {vi - Wi \ D~^Vi - Dl^Wi)g^ 



i=l 



> ^i\\Cx - Cyf + J]M,a;i||Z)ri^;, - Dr^Wifg^ 

i=l 

> (3(\\Cx - Cyf + f2u:,\\D-\ - DT^w.Wl) 

^ i=i ^ 

= /3||Q(x,Ul, . ..,Vm)- Q{y,Wi,.. .,Wm)\\K- 



(3.12) 



Therefore, by (1.1), Q is /3-cocoercive. It is shown in [16, Eq. (3.12)] that under the condition (3.1), 
zer(M + 5 + Q) / 0. Moreover, [16, Eq. (3.21)] and [16, Eq. (3.22)] yield 



(x, v) £ zer(M + S + Q)^xe'P and v gV. 



Now, define 



{x,vi, ...,Vm)^ ir'^x - ^0JiL*Vi,a^^vi - Lix, . . . ,0-;;^?;^ - Lmxj. 

^ i=i ' 

Then V is self-adjoint. Let us check that V is p-strongly positive. To this end, define 

T-.U^ Q: x^ {^^yfalLix, ... , ^/omLmx] 

Then, 



(Vx G v.) \\Tx\\g = '^u)iai\\Lix\\g. < \\x\\'^'^ojiai 



I II -^J II 5 



i=l 



i=l 



which implies that 



Now set 



< Y,^i<^i\\Li 



1=1 



.1 



m 

T ^fJjWiii^ji 
i=l 



Then, it follows from (3.2) that 5 > 0. Moreover, (3.17) and (3.18) yield 

m 

T\\Tf{l + 6)< r(l + 6)Y,^^^^\\Li\\^ = (1 + ^)"'- 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



(3.19) 



i=l 
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For every x = {x,vi, . . . , Vm) in K, by using (3.19), we obtain 

m m 



T-^lk||2 



r-i||x||2 



i=l i=l 
m m 

i=l i=l 
m 

2 



ix I ^/a7 



j=i 

m 

1 11 ..n V ^ 1 II ||0 

> IfII +z^^i WiPilb, 



WTxWl. 



i=l 



> (l - (1 + 5)-^) uAn{T-\a^\ ...,a'^ 



Therefore, V is p-strongly positive. Furthermore, it fohows from (3.20) that 

V'^ exists and < p'^. 

(i): We first observe that (3.3) is equivalent to 

T~^iXn - Pn) - Eitl ^iL*Vi^n " CXn G 

-z + A{pn - a2,„) + ai,„ - r-^a2,n 

for i = 1, . . . , m 

<^l^{Vi,n - qi,n) - Li{Xn - Pn) " D^^Vi^n G 



(Vn G N) 



Now set 



(Vn G N) 





— (Xj2, fl^ri) ■ ■ ■ 




Vn 


= {Pn,qi,n,--- 


) Qm,n) 




= {o-2,n, bl^n, ■ ■ 


■ 5 bra,n 




— (o-l^ni Ci^ni • • 


• ) Cm,n 




= (T~"^a2,n,C^r 





We have 



||an|k < +00, ^ llcnllx; < +00, and ^ ||dn||/c < +oo. 

nSN nSN nSN 

Furthermore, (3.22) yields 

(Vn G N) ~ ~ '^^^ G (M + S')(y„ - a„) + ^a^ + c„ - d„ 



Next, we set 

(VnGN) bn = V-\{S + V)ar, + Cn-dn). (3.26) 

Then (3.24) implies that 

WKWk < +00. (3.27) 

Moreover, using (3.21) and (3.26), we have 

(Vn G N) V{xn - Vn) - Qxn S (M + 5)(y„ - a„) + 5a„ + c„ - d„ 
4^ (Vn G N) {V- Q)xn G (M + 5 + F)(y„ - a„) + (5 + F)a„ + c„ - d„ 

4^(VnGN) y„= (M + 5 + F)~'((F-Q)a;„-(S + y)a„-c, + d„) + a„ 
4^(VnGN) y„= (ld + F-i(M + 5))"^((ld-F-ig)a;„-b„) +a„. (3.28) 
We derive from (3.25) that 

(VnEN) Xn+i = Xn + \n{{ld + V-^{M + S)y^{xn-V-^Qxn-bn) + a„ - x^) 

= Xn + \n (yJA {Xn " BXn - bn) + a„ - Xn^ , (3.29) 

where 

A = V-^{M + S) and B = V-^Q. (3.30) 

Algorithm (3.29) has the structure of the forward-backward splitting algorithm [13]. Hence, it 
is sufficient to check the convergence conditions of the forward-backward splitting algorithm [13, 
Corollary 6.5] to prove our claims. To this end, let us introduce the real Hilbert space JCy with 
scalar product and norm defined by 



{y{x,y)eKxK) {x\y)y = {x\Vy)^ and \\x\\v = ^J {x \ Vx),^, (3.31) 
respectively. Since V is a bounded linear operator, it follows from (3.24) and (3.27) that 

^ ^ llctnllv < +00 and ^ ^ W^n \\v < +00. (3.32) 

Moreover, since M -|- is monotone on JC, we have 

(V(a;, y) e K X K) {x - y \ Ax - Ay)^ = {x - y \ V Ax - VAy)-^ 

= {x-y\{M + S)x-{M + S)y)^ (3.33) 
> 0. (3.34) 

Hence, A is monotone on ICv- Likewise, B is monotone on /Cv- Since V is strongly positive, 
and since M -|- S is maximally monotone on /C, A is maximally monotone on K,v- Next, let us 
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show that B is (/3/>)-cocoercive on Kv- Using (3.12), (3.20) and (3.21), we have 



(V(a;, y)eKvX Kv) {x - y \ Bx - By)y = {x - y \ VBx - VBy)^ 

= {x-y\Qx- Qy)j^ 



>/3| 
= /3| 
= /3| 
>P\ 
>/3| 
= P\ 

= /3| 
> 



Qx - Qy\\K 

Qx - Qy\\K\\Qx - Qy\\K 



-1 



V '-WWQx - Qy\\jc\\Qx - Qy\\jc 



V 'Qx - V-'Qy\\K\\Qx - Qy\\K 
{V-^Qx - V^Qy I Qx - Qy)^ 
1 {Bx -By\Qx- Qy)^ 



Bx-By\\\r 



Bx - By\ 



(3.35) 



Hence, by (1.1), B is (/3/9)-cocoercive on K-v- Moreover, it follows from our assumption that 
2/3/3 > 1. Altogether, by [13, Corollary 6.5] the sequence {xn)n€N converges weakly in ICy to some 
X = {x,vi, . . . ,Vm) G zer(A + B) = zev{M + S + Q). Since V is self-adjoint and V^^ exists, the 
weak convergence of the sequence {xn)n&i to x in K-v is equivalent to the weak convergence of 
{xn)n£n to X in /C. Hence, x^ ^ x & zer{M + S + Q). It follows from (3.13) that x £ V and 
(jJi, . . . ,Vm) G T^- This proves (i). 



(ii)&:(iii): It follows from [13, Remark 3.4] that 



nGN 



Bx\ 



V 



< +00. 



On the other hand, from (3.20) and (3.36) yield Bxn — Bx = V ^{Qxn — Qx) 
that Qxn — Qx —7- 0. Hence, 

Cxn —J" Cx and (Vi e {1, . . . , m}) D^^Vi^n D~^Vi. 



(3.36) 
0, which implies 
(3.37) 



If C is uniformly monotone at x, then there exists an increasing function </>c : [0, +00 [ — )■ [0, +00] 
vanishing only at such that 



4'ci\\Xn - x\\) < {Xn- x \ CXn - Cx) < \\Xn - x\\ \\CXn - Cx\\. 



(3.38) 



Notice that (x„ — x)„gN is bounded. It follows from (3.37) and (3.38) that This proves (ii), 

and (iii) is proved in a similar fashion. □ 

Remark 3.2 Here are some remarks concerning the connections between our framework and 
existing work. 



(i) The strategy used in the proof of Theorem 3.1(i) is to reformulate algorithm (3.3) as a 
forward-backward splitting algorithm in a real Hilbert space endowed with a suitable norm. 
This renorming technique was used in [22] for a minimization problem in finite-dimensional 
spaces. The same technique is also used in the primal-dual minimization problem of [18]. 
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(ii) Consider the special case when z = 0, and {Bi)Ki<m and (-Dj)i<i<m are as in (1.4). Then 
algorithm (3.3) reduces to 

(Vn G N) Xn+l = Xn + Xn (^JtA (^Xn - T{CXn + ai,n)) + a2,n - , (3.39) 

which is the standard forward-backward splitting algorithm [13, Algorithm 6.4] where the 
sequence (7n)neN in [13, Eq. (6.3)] is constant. 

(iii) The inclusions (1.7) and (1.8) in Example 1.3 can be solved by [7, Theorem 3.8]. However, the 
algorithm resulting from (3.3) in this special case is different from that of [7, Theorem 3.8]. 

(iv) In Problem 1.1, since C and (-Dj"^)i<i<m are cocoercive, they are Lipschitzian. Hence, 
Problem 1.1 can be solved by the algorithm proposed in [16, Theorem 3.1], which has a 
different structure from the present algorithm. 

(v) Consider the special case when z = and (Vi G {1, . . . , m}) Gi = Ti, Li = Id, = 0, = 0. 
Then the primal inclusion (1.2) reduces to 

m 

find X such that £ Ax + '^ i^iBiX + Cx. (3.40) 

i=l 

This inclusion can be solved by the algorithm proposed in [26], which is not designed as a 
primal-dual scheme. 



4 Application to minimization problems 

We provide an application of the algorithm (3.3) to minimization problems, by revisiting [16, 
Problem 4.1]. 

Problem 4.1 Let T-Lhe a. real Hilbert space, let z E let m be a strictly positive integer, let 
{uji)i<_i<rn be real numbers in ]0,1] such that Yll^i^i — ^i f ^ ^oi'H)^ and let h: 7i ^ W 
be convex and differentiable with a /i~ ^-Lipschitzian gradient for some // G ]0,+oo[. For every 
i £ {1, . . . ,m}, let Gi be a real Hilbert space, let G Qi, let gi E ro(^i), let ii £ To^Qi) be 
fj-strongly convex, for some E ]0,+cxd[, and suppose that Li: % ^ Qi \s a, nonzero bounded 
linear operator. Consider the primal problem 

m 

minimize f{x) + > uJi{gi □ £i){LiX — Vi) + h{x) — {x \ z) , (4.1) 

i=l 

and the dual problem 

z -Y^uJiL*Vi] +Y^uji{g*{vi)+ti{vi) + {vi \ ri)g). (4.2) 
We denote by V\ and Pi the sets of solutions to (4.1) and (4.2), respectively. 
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Corollary 4.2 In Problem 4^.1, suppose that 



E ran Uf + □ dli){Li ■ -n)) + Vh\ . 

^ i=i ' 



(4.3) 



Let T and {(Ti)i<i<m be strictly positive numbers such that 



1 - T^^aiUJi\\L. 
\ i=i 



(4.4) 



Let e G ]0, 1[ and let {\n)nen be a sequence in [e,l\, let xq G H, let (ai^„)„gN cind (a2,n)neN be 
absolutely summahle sequences in %. For every i G {1, . . . ,m}, let Vifi £ Qi, and let {bi^n)nen (md 
(cj,n)neN be absolutely summable sequences in Qi. Let {xri)ne'H and (fi,n> • • • ,Vm,n)nm be sequences 
generated by the following routine 



(Vn e N) 



Pn = prox^j \^Xn - EilLl ^iL*'"i,n + V/i(j;„) + ai^n ~ ^)) + "-V 

Vn — Xfi 

— ~l~ ^nijPn ^n) 

for i = 1, . . . ,m 

Then the following hold for some x gVi and (vi, . . . ,Vm) G 'Di- 



(4.5) 



(i) Xn^ X and {vi^n, ■ ■ ■ ,Vm,n) ^ (^-l, • • ■,Vm)- 

(ii) Suppose that h is uniformly convex at x. Then 

(iii) Suppose that £* is uniformly convex at vj for some j G {1, . . . , m}. Then Vj^. 



Vn. 



Proof. The connection between Problem 4.1 and Problem 1.1 is established in the proof of [16, 
Theorem 4.2]. Since Vh is //"^-Lipschitz continuous, by the Baillon-Haddad Theorem [4, 5], it is 
/i-cocoercive. Moreover since, for every i G {1, . . . ,m}, is i/j-strongly convex, dii is i^j-strongly 
monotone. Hence, by applying Theorem 3.1 (i) with A = df, JrA = prox^j, C = V/i and for every 
i G {1, . . . D^^ = V£*, Bi = dgi, J^.^-i = prox^.^. , we obtain that the sequence (xn)nGN 
converges weakly to some x £ T-L such that 



G df(x) + J2 ^iL* {{dgi □ d£^){LiX - r,)) + V/i(a 



(4.6) 



i=l 



and the sequence {{vi, 

,ni ■ ■ ■ 1 ^m,n))ngN Converges weakly to some (f i, . . . , Vm) such that 

- TZl ^^^^^ e df{x) + V/l(x) 



{^xeU) 



(Vi G {1, . . . ,m}) Vi G {dgi □ dli){LiX - r^). 



(4.7) 



As shown in the proof of [16, Theorem 4.2], x G T^i and (iJi,...,!;^) £ ^i- This proves (i). 
Now, if h is uniformly convex at x, then V/i is uniformly monotone at x. Hence, (ii) follows from 
Theorem 3.1(ii). Similarly, (iii) follows from Theorem 3.1 (iii). □ 
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Remark 4.3 Here are some observations on the above results. 



(i) If a function ip: 7^ — )• M is convex and differentiable function with a /3~^-Lipschitzian gradi- 
ent, then V93 is /3-cocoercive [4, 5]. Hence, in the context of convex minimization problems, 
the restriction of cocoercivity made in Problem 1.1 with respect to the problem considered 
in [16] disappears. Yet, the algorithm we obtain is quite different from that proposed in [16, 
Theorem 4.2]. 

(ii) Sufficient conditions which ensure that (4.3) is satisfied are provided in [16, Proposition 4.3]. 
For instance, if (4.1) has at least one solution, and if ?^ and {fdi)\<i<m are finite-dimensional, 
and there exists a; G ri dom / such that 

(Vi € {1, . . . , m}) LjX — Tj G ridom^ij + ridom^j, (4.8) 

then (4.3) holds. 

(iii) Consider the special case when z = and, for every i G {1, . . . , m}, = 0, o"j = cr G ]0, +oo[, 
and 

U:v^\ : (4.9) 
1+00 otherwise. 

Then, (4.5) reduces to 

Pn = prOX^j (^X„ - T Xllll ^iL*Vi^n + ^KXri) + ai,n) ) + a2,n 
Vn — 2pji Xn 
Xn+1 — Xn -\- XniPn Xn) 

for z = 1, . . . , m 

qi^n = prox^g. [vi^n + cr{Liyn + Ci,„)^ + 

which is the method proposed in [18, Eq. (36)]. However, in this setting, the conditions (4.4) 
and (4.3) are different from the conditions [18, Eq. (38)] and [18, Eq. (39)], respectively. 
Moreover, the present paper provides the strong convergence conditions. 

(iv) In finite-dimensional spaces, with exact implementation of the operators, and with the fur- 
ther restriction that m = 1, /i: x 1— t- 0, £1 is as in (4.9), ri = 0, and z = 0, (4.5) remains 
convergent if A„ = A G ]0,2[ under the same condition presented here [22, Remark 5.4]. If 
we further impose the restriction = 1, then (4.5) reduces to the method proposed in [10, 
Algorithm 1]. An alternative primal-dual algorithm for this problem is proposed in [12]. 

Acknowledgement. I thank Professor Patrick L. Combettes for bringing this problem to my 
attention and for helpful discussions. 



(Vn G N) 
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